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Abstract. The main result in this paper is a one term Szego type asymptotic 
formula with a sharp remainder estimate for a class of integral operators of 
the pseudodifferential type with symbols which are allowed to be non-smooth 
or discontinuous in both position and momentum. The simplest example of 
such symbol is the product of the characteristic functions of two compact sets, 
one in real space and the other in momentum space. The results of this paper 
are used in a study of the violation of the area entropy law for free fermions 
in 1181 . This work also provides evidence towards a conjecture due to Harold 
Widom. 



1. Introduction 

The problem of estimating entanglement entropy (EE) is currently of consider- 
able interest in the physics community, in particular in condensed matter physics 
and in the theory of quantum information. The interest in EE in condensed mat- 
ter systems is due, in particular, to its scaling behavior and universal properties 
near quantum phase transitions, see for example 1251 \'M\ I!i2| and in particu- 
lar [AH\ . EE is an accepted measure of entanglement: in the quantum informa- 
tion context, entanglement is necessary for performing quantum computations, see 
e.g. |51 1441 I3l)[ |77) . An experimental demonstration of entanglement effects 

in a macroscopic system was reported in |14| . A connection between EE for spin 
chain models and Random Matrix Theory was found in [231 . 

The evaluation of EE in physical systems of interest presents considerable math- 
ematical difficulties ^1 153 El ^] ■ In several models EE turns out to be closely 
related to various versions of the strong Szego limit theorem (SSLT) for dimension 
one [52] and also for the higher dimensional case ^SIj see Remark |2 . 41 below . More 
precisely, the asymptotic behavior of EE as the size of the subsystem of interest 
becomes large for the XX spin chain model with a transverse magnetic field was an- 
alyzed rigorously in 22 using a certain Fisher-Hartwig theorem for large Toeplitz 
determinants established in , see also [2 (self-correlations for a translation 
invariant spin chain can be expressed in terms of Toeplitz determinants). For the 
more general XY model with a transverse magnetic field, the authors in |21| used 
the Riemann-Hilbcrt approach and the steepest descent method to find an explicit 
expression for EE (in the asymptotic regime where XY XX, the results in [23 
are recovered from the corresponding expression in j21|). 
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Note that the above mentioned results are for the one dimensional case only (we 
refer the reader e.g. to the references in |21l KMj for further results concerning spin 
chains and also harmonic chain systems). Much less is known about the higher 
dimensional case. It was shown in |34| that for the harmonic lattice model the 
entaglement entropy of a cubic region of size A behaves like A''^^, A — > cx), where 
d > 1 is the dimension, i.e. the EE is of the order of the area of the boundary 
of the cube. In the physics community this type of behavior is referred to as the 
area law for the entropy. The area law was initially discovered in the context of 
the so-called geometric entropy, which is a component of the Bekenstein-Hawking 
black hole entropy, sec 6 , 35 and the references in Qi^ . 

It is known that there is a correspondence between ID spin models and sys- 
tems of non-interacting fermions on a ID lattice by means of the Wigner- Jordan 
transformation. It is therefore natural to consider EE for fermionic systems in 
higher dimensions. It was noticed recently |18[ I43j . that the EE for a system of free 
fermions of arbitrary dimension on a lattice or in the continuum at zero tempera- 
ture violates the area law and is, in particular, of larger order than the boundary 
area of the region in which the entropy is evaluated. In |18j the authors make a 
connection between the asymptotics of EE and a conjecture due to Widom <4HI42| 
(see (12 .141) below), and then utilize this connection to posit a formula for EE of a 
continuous system, see H1.4f) below. The conjecture of Widom appeared originally 
in the context of time-frequency limiting problems, i.e. problems that involve the 
extraction of information about a signal from a measurement in a finite time — finite 
frequency window. 

In USIj some of the results are presented without proof. In this paper, we study 
Szego type asymptotics for operators of pseudodifferential type with discontinuous 
symbols in the higher dimensional case: Various specializations of these asymptotics 
provide the proofs of most, but not all, of the results left open in JS|i as explained 
below. 

To fix ideas we recall first some basic facts concerning EE (see e.g. |5j). Let 
Ha, Hb be two Hilbert spaces, which we assume for simplicity to be finite dimen- 
sional. Using the Schmidt decomposition 33 , any state tp S Ha <Xi Hb can be ex- 
pressed in the form t/j ~ J^i Ci(t>A.i®<t>B,i, where < q < 1, cf — I and (f>A,i, (t>B,i 
are orthonormal in Ha , Hb , respectively. Associate with ijj the density matrix 
p = ^ip* acting on Ha ® Hb- The reduced density matrix pA acting on Ha is de- 
fined as PA = Tr/fa p^J2i Ci(f>A.i(f>*A.i and similarly pB = Trn^ P^J2i (^i^B,i(f>*B,t- 
It is easy to check that the matrices pA and pb are well-defined independent of the 
choice of the orthonormal vectors (t>A,ii 4'B.j- 

The entanglement entropy of a state ip € Ha ® Hb measures how far the state 
■0 is from a product state of the form (pA ® 4>b^ and is defined as the von Neumann 
entropy of either of the reduced density matrices 

S = - Tt{pa log2 Pa) = - Tr(pi3 log2 Pb) = - ^cj log2 cj 

i 

which is precisely the Shannon entropy of the (squared) Schmidt coefficients 

In many problems, one is interested in finding EE for tp which is the ground state 
of some general many body system. In |26| the EE of a system of non-interacting 
(free) fermions in the ground state was studied. We note that the case of a general 
system with interactions is very difficult and at present time very little seems to be 
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known for dimensions higher than one (see however !13| and !34' where the area law 
is derived for the Kitaev model and for the harmonic lattice model, respectively). 

Let r C M'' (resp., F C T'') denote a compact set in momentum space and fix a 
compact region 57 C M'^ (resp., C Z'') in position space for the continious (resp., 
lattice) case. The ground state of the system is defined by the projection P in 
L'^{W^) (resp., P{'L'^)) onto the modes in the Fermi sea F. We study the entropy 
of entanglement between fermions located in a compact region 17, scaled by some 
large A, in position space, and its complement. Let Q — x\n be the projection onto 
\n in L2(Krf) (resp., 

For the system at hand, all the important quantities for the entanglement prob- 
lem can be described in terms of the operator PQP |2f)| : in particular the average 
number of fermions in Af2 is given by <N> = Tr PQP, the particle number variance 
is 

(1.1) (AiV)2 = Tr[PQP(l - PQP)] 
and the EE is given by 

(1.2) S = Sn,rW = Trh{PQP) 
where 

(1.3) h{t) = -t l0g2 t-{l-t) log2(l - t). 

Note that P = Txv^'^^ where T denotes either the Fourier transform or the 
Fourier series in the continuous, lattice case, respectively. 

Assuming the applicability of the Widom conjecture (|2.14() below to the function 
(|1.3() the authors in suggest the following explicit leading order asymptotics for 
the EE of a continuous system as A — > oo 

(1.4) ^n,r(A)= 'IS^ l^ / / \n.-nMS.dS^ + o{\''-^\og^\\ 

(27r)" 12 Jq^ Jqy 

where rix^n^ are outward unit normals to the (smooth) boundaries 951, 9F and 
dSxi dS^ are the area elements. Although the formula (|1.4|) was conjectured in [TH] 
only for continuous models with smooth boundaries, it is probably also true for 
piecewise smooth boundaries in both the continuous and the lattice case. Indeed, 
(|1.4I) was recently checked numerically in the lattice case for d = 2 in [3j and for 
d = 2, 3 in |28| . and an extremely close agreement concerning both the order and 
the (leading) coefficient was found. 

The following two results proven in |18) provide corroborating evidence towards 
(jT^ . For the cubic domains F = [-i, i]'' and n = [0, 1]'' (resp., n = {0, 1}'*) in 
the continuous (resp., lattice) case, the following holds as A ^ oo 

(1.5) l(A)"V(A)<^.,r(A)<d(Ay"^5,(A) 

where S'i(A) is the entanglement entropy for the one dimensional system with Fi = 
[— i, i] and ill = [0, 1] (resp., {0, 1}). We note that H1.5|l together with a result in 
[221 for the lattice case 

(1.6) 5i(A) = ilog2A + o(log2A), A^oo, 
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implies 

(1.7) K^)'"^°^^'-^^^^^'^-K^)'"^°^^' 

thereby demonstrating the violation of the entropy area law for cubic domains in 
the lattice case in a way which is consistent with (|1.4|l . 

In fTSI the authors also prove that for arbitrary (measurable) compact il,r in 
both the continuous and lattice cases 

(1.8) 4(AiV)2 < 5oj(A) < C ■ (log2 A) • (AN)^ 

where the constant C depends only on the dimension d. (The estimate H1.8|l for 
the lattice case with d = 1 was proved in ^^l- The proof of H1.8|l in ^| for all 
d > 1 in the lattice case is analogous, but the continuous case requires a new idea 
as provided in J 8;-) 

The fact that Theorem 12.11 below is sharp implies, together with ^1.1^ . that in 
the continuous case if F, ft have boundaries then for some ci, C2 > 

(1.9) ciA'^-MogaA < (AiV)2 < caA'^-MogaA, A ^ oo, 

and also that for any /3 e (0, 1) there exists (a Cantor-like) set F such that for 
ri = [0, 1]'' for some ci, C2 > 

(1.10) ciX'^-'^ <iANy <C2X'^~^, A^oo. 

The inequalities (|1.8|) . (|1.9|l yield the following: In the continuous case, if il, F have 
boundaries then for some ci, C2 > that depend on 17, F 

(1.11) ciA'^-Mog2A < ^o,r(A) < C2A'^-i(log2A)2, A ^ oo. 

Note that Hl.ll|l proves the violation of the area law in the continuous case for 
arbitrary domains with boundary and also gives the expected order for the 
lower estimate (and also the expected upper estimate, up to a power of log2 A) 
consistent with lll.4|l . 

Remark 1.1. It was proved independently in '43' that for cubic domains in the 
lattice case the estimate Hl.ll|l holds. Note that in this case we have the stronger 
estimate (jl.7|l . 

Concerning irregular boundaries, H1.8|l together with (|1.10() shows that in the 
continuous case in any dimension d and for any (3 G (0, 1), there exists a set F so 
that for n ~ [0, 1]'' for some Ci, C2 > 

(1.12) ciX'^-l' < 5n,r(A) < caA'*-'' log2 A, A ^ oo. 

Finally, we note that in the one dimensional lattice case the following result is proved 
in For any /3i G (0, 1) there exists a set Fi C T such that for fli = {0, 1} for 
some ci , C2 > 

ciX^-^' < Sn,,rA^) < C2X^-^' log^ A, A ^ oo. 

Combining this result with a straightforward generalization of (|1.5(l to the case 
T = Tf, fl — i^f, we see that ()1.12f) holds also for the lattice case (where /? = l—fii)^. 

^The fractal set F in [HI is very similar to the fractal set appearing in Lemma l2.4l below. which 
was constructed in earlier work of the author, see math. FA/0212215, math.CA/0212254. 
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The outline of the paper is as foUows: In Section 2, we describe the main results, 
the proofs of these results together with various technical lemmas are given in 
Section 3. 

2. Main results 

We now state the main results of the paper. As noted above, some of these 
results were used in ^2] to estimate EE for certain physical systems. 

The problem of evaluating the Szego type asymptotics for operators of pseudo- 
differential type with symbols discontinuous in both position and momentum in the 
higher dimensional case was introduced in 02 (see Remark [2.31 below) . we refer 
also to 1301 12ZI where certain related results can be found. See e.g. for a more 
detailed account of the results in this area. 

Let d gN, let mes denote the Lebesgue measure in 
II • life be the standard norm in L'^(R^), fc = 1, 2, oo. Denote by || • ||, || • Uss: II ' llsi the 
operator norm, the Hilbert-Schmidt and the trace-class norm in if, respectively. 
Let 5(R^) stand for the Schwartz space and denote by J the integration over M . 
Let and " denote the Fourier transform: !Fx->u[g{x)] = g{u) — J e^^"^'^ g{x) dx, 
g G H. For a function g G denote its modulus of continuity by uj2[g]{h) = 

||f;(- + h) — g{-)\\2, h G W'-. For a set C W- denote by xo. its characteristic 
function. We characterize the regularity of the set in terms of a-'2 [xn] : assume 
that xn S L^(M'^) (i.e. mes(r2) < oo) and that there exist < /3o < 1 and co > 
such that for small enough \h\ 

(2.1) {oJ2[xn]{h)f <cn\hf" 

where \h\ = {Y.%i Let Vth = [x - h\x e 17}, e R^'. The left-hand side 

of H2.1|l equals 

(2.2) ||xn(- + h)- xn{-)\\l = mcs{n \ Q,,) + mes(r! \ Q^h), 

which gives the geometrical meaning of (|2.1|) . Introduce a projection P : H ^ H 
by {Pg){x) = xn{x)g{x), x G K*^, g € H. Following 03 we consider a family of 
integral operators A\ : H ^ H , X > 2, oi pseudodifferential type with a non-smooth 
or discontinuous symbol. Let Ax, A > 2, have the kernel 

(2.3) Ka, (x, 2/) = (^) / e'^^-^^-y^a{x, y, 

where a (is measurable and) satisfies the following mild condition: Define 

(2.4) (/)(u) = sup |j^j^„[cr(x,y,C)]|, ip{u) = <j){u)(t){-u) 

x:,y 

and assume ip G L^(M'') and for certain < /3 < 1 and c > 

(2.5) / ■>p{u)du <cp-^, p>l. 

J\u\>p 

See Remark l2 . SI below for an example of such a. Setting A = /i^^, /i > 0, in H2.3|l we 
can obtain semiclassical type asymptotics. We study the asymptotics of the trace 
of 

(2.6) Pf{PAxP)P - PI{Ax)P, 

as A ^ oo, for suitable functions /. For the particular choice a{x,y,S,) = Xr(C)7 
r C M'*, such a study is motivated by the following question: what can be said 
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about a function if its restriction to Q and the restriction of its Fourier transform 
to r are known, see 02! • An order sharp estimate for the trace of H2.6(l is found in 
Theorem 12.21 below in two settings; we assume either that / is analytic, or A\ is 
self-adjoint and / has a bounded second derivative. Theorem l2.2l is a generalization 
of the classical Szego limit theorem see e.g., [321 EZ| for a review of related 
results. Let log denote the natural logarithm. The following result is basic for the 
proof of Theorem 12. 21 

Theorem 2.1. Let j3, a, A\, (3n, fl be as above. Then there exist two constants 
c(ri, a) and C{Q, a) such that one has, for A > 2, 

(2.7) \\PA^\\l^<c{n,a)-\\ 
and 



(2.8) \\PA^{I-P)re^<C{n,a) 



X^-HogX, P = f3. 



and the same estimate holds for \\{I — P)AxP\\q^. The estimate H2.8|l is sharp 
on the described class of a, that is for a certain a (which can be chosen so that 
the corresponding A\ is self-adjoint), the reverse inequality to H2.8|l holds for some 
(different) constant C(il,,a). 

For an analytic on some disc f{z) = J2m=o ^mZ™ set 

oo 

(2.9) /,(z)= ^m(m-l)|c™|z'"-2 
and 

(2.10) 5(a) = sup 

The condition S{a) < oo holds for a wide class of cr, see Remark 12.51 below . 

Theorem 2.2. Let the assumptions of Theorem \2.1\ hold. Assume further that 
either 

(i) S{a) < oo and f is analytic on a neighborhood of {z : \z\ < S{a)}, or 

(ii) Ax is self-adjoint for all X > 2, and f is such that f" £ i°°(R). 

Then the operator H2.6|l is trace-class. Moreover, there exists a constant C{fl,a) 
.such that one has, for X > 2, 



(2.11) Tv[PfiPAxP)P-PfiAx)P] <Cin,a)-C{f)- 
where 



A'^-'^logA, P = pn, 



2 [||/"||oo, in case (ii). 

In both cases, the estimate H2.11(l is sharp, that is for certain f and a, the reverse 
inequality to (|2.11(l holds for some (different) constant C(Q.,cr). 

In the case of analytic / we prove slightly more, namely that the sharp estimate 
(|2.11|l holds with the trace-class norm of (|2.6(l in the left-hand side. See subsections 
13.11 and 13.21 for the proofs of Theorem 12.11 and Theorem 12.21 respectively. It is 
possible to compute the leading term in the asymptotics of Tt Pf{PA\P)P in 
certain special cases. 
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Corollary 2.3. (i) Let and T be two bounded domains with C°° boundaries. 
Assume that the symbol a does not depend on y, and that a[x,S^) = ''"(a;, ^)xr(C); 
where r G iS(M^''). Let f be analytic on a neighborhood of the disc {z : \z\ < ||t||oo} 
and satisfy /(O) — 0. Then the operator f{PA\P) is trace-class, and furthermore, 
for any small enough e > there exist two constants C'{fl,T,T, e) and A(e) > 2 so 
that, for A > A(e), 



(2.12) 



TvfiPAxP)- 



f{T{x,0)dxd^ 



r Jn 



< C{n, r, T, e) • ((1 + e)||T||oo) • A-^-i log A. 



(ii) Let (|2.1|l hold and assume that the symbol a is real-valued, depends only on 
^ and (ESI holds. Let f" G L°°([- ||cr||oo, iVlU]) and assume /(cr(0) £ L^iW^). 
Then f{PA\P) is trace-class and there exists a constant C(r2, a) such that, for 
A > 2, 



(2.13) 



Tr{Pf{PA^P)P)~(^^ mes{n) j f{a{0)d^ 

<c{n,a)-\\n\ 



The proof of part (i) follows from the functional calculus results developed in 021 1 
see subsection 1^21 Part (ii) is proved as follows. Under the above assumptions Ax 
is self-adjoint, A > 2, and also the operator f{Ax) is well-defined and has the kernel 
K^(^,)(x,y) = Je^xi-i-y) f{a{0)d^. The operator PfiAx)P, A > 2, is 

trace-class since it is a composition of two Hilbert-Schmidt operators PJ-~^\f\^^^J-, 
J^~^\f\^^^{sgn f)TP. Note also that K f(^ji^^{x,y) is continuous, since /(o'(^)) G 
L^(R'^). Now we simply write the trace of Pf{A\)P as the integral of its kernel 
over the diagonal. 

Remark 2.1. Let xsi,Xr satisfy (|2.16() below with < /3si,/3 < 1, respectively, and 
set cy{x, y,£) = xr(C) (this example is considered in the proof of Theorem l2.1|l . The 
order sharp remainder estimate in H2.13|) for ^ j3 shows that the set with less 
regular boundary contributes to the order of the remainder. In the case Pn = f3 
the logarithmic factor persists even for f2, F with C°° boundaries (in which case 
— (3 — I) due to the fact that the symbol Xfi{x) Xr(C) Xfiiu) has discontinuities 
in both the position variables x,y and the momentum variable ^. (This should be 
compared with the power type asymptotics in e.g. |27[I2L)|[TC| when a discontinuity 
in momentum only is present.) 

Remark 2.2. The leading term Tr Pf{Ax)P in the asymptotics of Tr P/(P^a-P)^' 
in (|2.12l) and (|2.13|l is of Weyl type. It is written as an integral over the diagonal 
(x,x,^) of the corresponding to A phase volume. 



Remark 2.3. In ^2 021 the following second order generalization of (|2.12|l was 
conjectured: Under the assumptions of Corollary 12. 3^ ) the following holds, as 
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o Jr 



\2ttJ 

(2.14) / A yi- HogX 



(t') / / \nx ■np\U{0,T{x,^);f)dSxdSi 

\2nJ 47r^ Jan Jar 

o(A''-MogA) 



where nx,n^ are the outward unit normals to 951, SF, respectively, and 

This conjecture is still open and it was one of the motivations for the present work. 

Remark 2.4. In a broader context, the Widom conjecture (|2.14(l is a generalization 
of the strong (two-term) Szego limit theorem (SSLT) for the continuous setting. (As 
noted earlier, a generalization of the SSLT plays a central role in the computation 
of the EE for the XX spin chain model in [22].) The SSLT was initially used by 
Onsager in his celebrated computation of the spontaneous magnetization for the 2D 
Ising model (see e.g. [3 and the references therein). It is interesting to note that in 
Onsager's computation (and also in \22\) the leading asymptotic term vanishes, and 
one needs to compute the sub-leading term. This is exactly the situation in [TH| : 
the leading term should vanish since for the entropy function h in /i(l) = 

(and also /i(0) = holds). 

Remark 2.5. Let F C M'^ be such that for some < /? < 1 the function xr satisfies 
(|2.1I) with J — (3 (see also Lemmas l2 .41 12 . 51 below) . Set a{x,y,£,) — T{x,y, ^) xriO^ 
where t is satisfies for some c = c(t) < oo 

(2.15) snp\T^^4T{x,y,C)]\<cil + \u\y''-\ ueR''. 

A standard application of the Cauchy inequality implies that there is C = C (r, d) < 
oo so that ipiu) <CJ \xriu-v)\'^{l + \v\)-''-^dv e L^{R'^), and (EHl holds (see, for 
instance, TH', Section 3.4.2]). For a{x,y,^) = xr(C) we can take V(w) = lxr(u)p. 
Now for H2.10|l . If a{x, y,£,) is a classical zeroth order (parameter dependent) symbol 
in the sense of pseudodifferential operators, then S{a) < oo is a standard result, see 
e.g. [inilig. However we are interested here in a with limited regularity for which 
no such general results are available. We restrict ourselves to the following standard 
example in which irregularities in x and, most importantly, in ^ are allowed. Assume 
a = (t{x,£) (we only need S[a) < oo when Ax is not assumed to be self-adjoint) 
is of the form t{x,£,) xiO- Then A\ = A\T~^x^ where A\ has integral kernel 
{^Y J e'^^'^-y'><T{x,^)d^. Assume that x e i°°(K'^), then \\T-^xH < llxlU 
(the factor x is allowed to be quite irregular). Now sup_)^>2 ||^a|| < oo follows 
from the standard estimate for the bilinear form under the sole assumption that 
(f){u) = snp^\J^^^uT{x,^)\ e L^{R'^). Indeed, by the Cauchy inequality for any 
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{2^Y\{AU,g)\<\'' jj ~^{-\{x-y))\f{y)\\g{x)\dxdy 

< [X'^jj 4^{-\{x-y))\f{y)\HxdyY^{\'' jj 4^{-X{x - y))\g{x)\Hxdy)"^ 

< Il0l|l||/I|2||<?||2 

which imphes S{a) < {2tt)^'^ I|0||i- Note that if cr{x,^) — T(a;,^)xr(0 where F is 
as before, and sup^ \ J^^^ut{x, ^)| < c (1 + as in H2.15|l . then both (|2.5|l and 

S{cr) < oo hold. 

In the proof of Theorem 12.11 we use the foUowing auxihary results. They are 
concerned with a two-sided version of the estimate H2.1|l : assume that g E L^(R'*) 
is such that there exist < 7 < 2 and ci, C2 > so that for small enough \h\ 

(2.16) ci|/iP < {u;2[g]{h))^ <C2\hp. 

Lemma 2.4. For any d G N and < /3si < 1 there exists a compact set fl C R''' 
such that xn satisfies (|2.1()|) with 7 = 

Lemma 2.5. Assume that a function f G L^(M'^) satisfies H2.16|l for some < 7 < 
2 and ci, C2 > 0. Then 

(i) there exist 61,62 > such that 

(2.17) b.p-'Kf \fm^d^<b2p-\ p>l; 

"'l?l>p 

(ii) the upper estimate in l|2.16|l implies the upper estimate in (|2.17|l . 

Lemma |2.4I is proved in subsection 13.41 Lemma |2.5l for 7=1 was proved in |H1 
Lemma 2.10, 4.2], the proof for 7 G (0,2) is analogous and is left to the reader, 
see also ^7] and ^| Lemma 3.4.1]. (If one introduces an average of u}2[g]{h) over 
II^IU ^ then the upper estimates in this modification of H2.16|l . and in (|2.17() 
become equivalent, and so do the lower ones, see [TniITT) .') 

3. Proofs 

3.1. Proof of Theorem 12. IL We denote A-independent constants (that may de- 
pend on ri, ct) by c^, fc 6 N. Let us consider ||PA>(/ — P)||62 only, the case of 
II (/ — P)A\P)\\(s^ is completely analogous. 

1. Let us prove first the upper estimate p.8|l . Using (|2.4() we obtain 
(3.1) 

^.2d 



|PAA(/-P)r6. = ( ^ ) / / / / e^^^<'-^-y^ e~^^^<^'y^ a{x,y,0<j{x,y,v) 



< 



X xn{x){l - xn{y))didydridx 

2d 



2ti 



^5--A(x-s) [cr(a;,y,^)cr(a;,y,7?)] xn{x) (1 - Xn(j/)) dxdy 

T?— *A{a; -y) 

V'(A(x - y)) xn{x) [1 - XQ.{y)] dxdy. 
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Changing variables x — y = 2a;', x + y = 2y' and dropping the primes we rewrite 
the right-hand side of H3.1|l in the form 

(3-2) (^^^ J ^{2Xx) dx J xo(x + y)[l- xni~x + y)] dy. 

Since xo takes on values or 1 only, we have 

Xn{x + y)[l-xn{-x + y)]dy^\\xn\\l~ / xn{x + y) xn{-x + y) dy 



= \ J[xn{x + y)^ xn{-x + y)f dy = ^ iuj2[xn]i2x)y 
Therefore (IXT)) and imply 

\\PA^{I - P)\\%^ < X^^ c, U{2\x) {u2[xn]{'^x)f dx 



(3.3) 

^X^^cJ [ +[ +[ )^p{2Xx){iU2[xn]{2x)fdx. 

\j\x\<l/\ Jl/\<\x\<l J\x\>l/ 

Using (|2.1|l and making a change of variables u ~ 2Xx we estimate the first integral 
in (|3.3|l as follows 

X^'^ci [ V(2Ax) {uj2[xn]{2x))^ dx < A^'^ca / V(2Aa;) \xf" dx 
<X'^-P''C3( iij[u)\uf'^ du<X'^-^<-^ CiUWi, X>2. 

J\ii\<2 

Next, noting that {uj2[xn]{'^x))'^ < 2\\xfi\\2 = 2mcs(fi), setting u — 2Xx and using 
(|2.5I) we estimate the third integral in 13.311 in the following way 

\2(i „ / „;. /ox^N /, . r — l^'o™^^2 ^ o „,„„/ON \2d 

(3.5) 



A^'^ci / V(2Ax) (w2[xo](2a;))2da; < 2mes(rj) A^'^ / i;{2Xx) dx 

x\>l J\x\>l 



<X'^C5 [ tl;{u)du<X'^-^ cq, X>2. 

J\u\>2\ 



Consider now the second integral in (|3.3|) . Set ^'(r) = J^d-i ^{r9) dSg, r > 1, and 
note that by (|2.5|l 

/•OC 

(3.6) / 'i/{s)s'^-^ds <cr-'^, r>l. 

J r 

Then the second integral in (|3.3|1 is estimated as follows 
(3.7) 

A^'' ci / V'(2Ax) (w2[xn](2x))2 dx < A''"''" cy / V(") Wf"" du 

/2A 
*(r) r'^''+'^-^ dr. 

Writing ^'(r) = — /^°° ^(s) s''"^ and integrating by parts we obtain 

/2A r /"OC 

^(^)j.Ai+d-i^,. ^^d-ft, 2'3" y ■^{s)s'^-^ds 

„ 2/30^^*^2 y yii{s)s'^'^ds + (3n J [J -^{s) s"^-^ dsj rf^"'^ dr 
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Discarding the negative term above and using (|3.6|l we obtain 



2A 



/ 1, Pn<f3 

I ^(^) ^fti+d-1 < ^d-/3n cg ■ 1 + < log A, /3o = P 

^ [aA^-^ /3o>/3 



< Cg 



A'^-'^logA, /3 = /3, 



Substituting the latter estimate in (|3.7|l and collecting the estimates (|3.4() . (|3.5ll we 
complete the proof of p. 8(1 . 

2. Let us now prove the sharpness of the estimate ((2.8|l . We choose arbitrary < 
Pn,f3< 1. By Lemma lT^ there exist two sets fl, F such that xn,Xr & L'^{M.'^) satisfy 
the two-sided estimate H2.16|l with the power f3n and f3, respectively. Set a{x, y, ^) = 
Xr(0 (then the corresponding Ax is self-adjoint) and let i^iu) = |xr(u)p. Apply 
Lemmal^fi) to xr- Then for some 61,62 > 0, bip^^ < J^^^^^ip{u) du < 62/3 '^ for 
p > 1. Let K be such that 62^"^ < ^ ^^^^ k > 1. Then 

(3.8) ^p-^< I ^iu)du, p>l. 

^ J p<\u\<K,p 

In place of (|3.1|l we now have an equality and therefore 



||P^a(/-P)|||, - A^^cio / V(2Aa;) {L,2[xn]i2x)f dx 

(3.9) 



>A2'^cio/ V'(2Aa:) (a>2[xo](2a:))2da; 

Jl/A<|x|<l 



Jl/\<\x\<l 

> X'^''^"ci2 I iP{u)\uf"du 

^2<|u|<2A 

where we have use the lower estimate in (I2.16|) for xn- Let A > k, set L = 
[log^ A] where [■] denotes the integer part of a number. Now we split the domain of 
integration as a union of concentric domains in the standard way (see e.g. [S]) 

A'^-*' C12 / xP{u)\uf'' du > A'*-*' C12 V / 'ilj{u)\uf'^ du 

•'2<\u\<2\ J2K'<|ti|<K-(2K') 

>A'^-*'ci2X^(2a^')^"--(2/«')"'' 

1=0 

L 



1=0 



where we have used ((3.8(1 . Considering the three cases when /3a is smaller than, 
equal to, and greater than f3 separately, and using the fact that k~^X < < X 
together with ((3.9(1 we conclude that for the operator A\ 



(3.10) I|PAa(/-P)|||, >ci4 



X^-^logX, p^Pn 
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for A > K. The sharpness of the esthnate ()2.8|l is proved. 

3. Finally we prove p. 7(1 . Analogously to ((3.1|l . ((3.2|l we obtain 

IIPAaIII, < (^^^ y" ^(2Ax)dx y" xn{x + v)xn{~x + y)dy 

< (27r)-2rf2-''||V'||inies(17) • A'', 
for all A > 2. This finishes the proof of Theorem l2.1l 

3.2. Proof of Theorem l2.2L It is important that in both cases the fact PA\ G &2 
implies that the operator 1(2. 6|l is trace-class, and also that in this case the absolute 
value of the trace of (|2.6|l can be estimated as follows: there exists a constant 
C(r2,cr) such that one has, for A > 2, 
(3.11) 

Tr [PfiPA^P)P-PfiAx)P] < Cin,cT)-Cif)-\\PAxiI~P)\\e, \\iI-P)AxP\\e,- 

In the case of self-adjoint A\, (|3.11|) was proved in |57| Theorem 1.2] (note that 
Uo<t<i Ua>2 spec A\ C M). In the case of analytic / the idea of the proof goes back 
to 00] • More precisely, denote Q = I — P and note that for any to e N, m > 2 we 
can write 

PA™P = PAx{P + Q)A';^-^P = P{AxP)A';^-^P + {PAx)QA';^-^P 
= P{AxP)Ax{P + Q)Al'-^P + {PAx)QA'l'-^P 

(3.12) = PiAxPfA'^-^P + {PAxfQA"^-^P + {P Ax)Q A^;^-^ P = • • ■ = 

m— 1 

= P{AxPT + {PAxr-'QA{P 

Also for j >2 

A{P - Ai-\P + Q)AxP - Ai-'PAxP + Ai'^QAxP 
= A{-\P + Q)AxPAxP + Ai-'QAxP 

(3.13) = Ai-^{PAxfP + A{-^Q{AxPf + A{''QAxP = ■■■ = 

= Ax{PAxy-'P + J2 AlQ{AxPy-''P 

k=l 

Substituting (|3.13|) in H3.12|l we obtain 
(3.14) 

m — 1 

P(Aa)"P - {PAxPr = J2 {PAxVQAlP 

m—1 m— 1 i"! 

= Y.^PAxr~'QAx{PAxy-^P+Y.^PAxr-'QY,AiQ{AxPy~''P 

j=l j=2 k=l 

There are terms on the right-hand side in H3.14|l each containing both 

PAxQ and QAxP. Hence 

(3.15) 

WiPAxPr ~ P{AxrP\\^^ < ""^"1" ^h \Axr~^\\PAxQ\\e. WQAxPWe,- 
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Recall (imi and note that (pnT|l (even with \\Pf{PAxP)P - Pf{Ax)P\\^^ in the 
left-hand side) holds for any f{z) analytic on a neighborhood of {z : \z\ < S{a)}. 
Note finally that for /(z) = and a self-adjoint Ax 

Tr[Pf{PAxP)P-PfiAx)P] = \\PAx{I ^ P)\\l,. 
This together with the sharpness in Theorem 12 . II implies sharpness in Theorem l2.2l 

3.3. Proof of Corollary 12.31 We need to prove the case (i). By Theorem 12.21 
we only have to analyze Ti{Pf{Ax)P), f analytic. We use the functional calculus 
developed in |42) . Denote by opr the depending on the parameter A > 2 operator 
with integral kernel {2n)-'^ J e'^-'-^'-y^ t{x,£,/ X) d£, and let Qxr = opxAr(0- Then 
Ax = (opr) Qxr- By 021 Lemma D(iii)] for any e > there is a constant A(e) > 2 
such that for all A > A(e) 

llA.||<(l + 6)||r|U. 
Next, by 02 Lemma 2], for any e > there exist C(e) and A(e) > 2 such that for 
all m G N and A > A(e) 

(3.16) \\P{AxrP-P{opT"^)QxrP\\e,<C{e){l + er^\\T\\Z-X''-'. 

Clearly Tr(P (opr") Qal -P) = {^)'^ J^. J^{t{x,0)"' dx d£,. By the argument 021 
p. 184], as A — > cxD 

Tr(P(/(opr))Q;,rP)= (^^) jJj{T{x,0)dxdi + O{X'-^). 

Now for any f{z) = Ylm=i analytic on a neighborhood of {z : \z\ < ||t||oo} 

we have c,„ = 0((1 + ^) ™||r||™) for some 6 > 0. Taking e < S in (jXTSIl finishes 
the proof of Corollarv l2.3l 

3.4. Proof of Lemma l2.4l We start with d = 1. Here we write [3 — (3n for brevity 
and denote the Lebesgue measure in R by mesi. The case /? = 1 is trivial, any finite 
interval with non-empty interior would do. Let us therefore assume < /? < 1. We 
construct a set 57 with the required properties as a (finite union of sets each of 
which is a) countable union of closed intervals obtained by a process similar to the 
construction of a Cantor set. The difference is that we do not remove but rather 
add intervals. Let < a < oo, = [0, a] and < q < 1. We explain how q is 
related to the given (3 later. We construct the set fta as follows. We start with 
an empty set and at the 0th step we add the middle qth part of la to fla- Each 
of the remaining intervals has length aQ, where we denote Q = Note that 
< Q < i. Then we take the middle qih part of each of the two remaining intervals 
and do not add it to the set ^l^- Each of the four remaining intervals has length 
aQ^ . Now we take the gth middle part of each of the four intervals and add it to 
fie,. This completes the 1st step of the construction of SIq,. We continue in this 
manner. The set Via we obtain has the following properties: 

1. For each fc G N, it contains 2^*^ intervals of length aqQ^''; 

2. Half of these 2^'° intervals will have an interval of length aqQ^''~^, at a 
distance of aQ^'^'^^ to the right, which does not contain points from ila- 

For any fc G N introduce notation 

4"^ = aiqQ^'^ + Q2^-+1) = aiq + Q)Q^'' 
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and consider the shift of the set fla to the right by 

(3.18) he[ai''\bi"'>], keN, 

units. RecaU that Q e (0, 1/2) and set /3 = 1 - j^^jq e (0, 1). (Actually we start 
with a given /3 e (0, 1) and after that define Q as above and q — 1 — 2Q e (0, 1).) 
By the properties 1 and 2 above for h as in (|3.18|l 

1 1 , _ log 2 

(3.19) mesi{in^~h)\n^)>--2'"'-aqQ^'' = -aq-{Q^''y 

We note that in H3.18|l . ciQ^*^ < h < C2Q^^ , where ci = a[q + Q) and C2 = 
a[qQ^^ + Q) do not depend on k. In particular Q^*^ > a{qQ-^+Q) ^ ^'^^ ^ 
(I3.18|l . and hence H3.19|) implies 

(3.20) --^(("-'^)\^")^^R^Q^?TW-'' 

for h as in (I23HI), and therefore for h G [JkefMk'^ ^ 4"^] • Using q = 1 - 2Q we find 
that 

% = ^^^ = Q-Q^e(o,i/4), ,E(o,i), 

and hence the set Ufeewl'^i^^ ^i"^] gaps. We need (|3.2()|l to hold for all small h. 
In order to fulfill this condition, we consider a finite union of the scaled sets Via for 
appropriate a. More precisely, the first set in the union is the set Qi corresponding 

to a = 1. Then holds for h G [jkeni'^k^ '^k^]- ^0*^ ^^^^ ^'^^io 

is independent of fc and a. Now for this 7 we construct the set fl^. For the set $7^, 
f^lUli holds (wfth a different constant) for h e UfceN^ai^^ 7^1^^] = Uketii^k^ ^^^k^] 
because 70^"'' ~ b\^'' for all fc e N and a > 0, in particular for a = 7. Now for the 
union fii U {il-y — 2), 13.20() holds (with a different constant) for 

^ ( U u ( U [^i^7&i^^]) - U [4^76W]- 

^fceN ^ ^feGN ^ feeN 

Now we construct the scaled sets 51^2, • ■ • , fl^N , where the (finite) number G N 
is determined by the (independent of k and a) condition 7^^'^"'' > '^i-'i or 7^ ^ 
-^^^ (recall that 7 > 1 by JS^TJl). Set finally 

n = niu{n^-2)u {n^2 - 3 - 7) u • ■ • u fn^N - ^(i + 7^'"^)") . 

Since is finite, H3.2U|) for the set holds (with a different constant) for h in 
U ^] U U [7fei^7^^«] U • • ■ U [7~-6i^7"5W]) 

00 
fe=2 
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Therefore (|3.2()|l holds for the set fl with some constant for all h G (0, a^^'']. 

Now we explain why the estimate opposite to H3.20|l holds for the same /?. Again 
since TV is finite it suffices to consider one set Via with an arbitrary a > 0. Choose 
first any /c = 2, 3, • ■ • and h G [al^\ '^fc-\]- We estimate the contribution of the sets 

up to generation k to mesi((r2Q — h) \ fi^) by a^"_\ times their total number. We 
estimate the conribution of the sets from generation k + 1 and onwards by their 
total length. As a result we obtain the following estimate 

k oo 

mesi {{n^ - h) \ n^) < ai"\ ^ 2^' + ^ aqQ^' ■ 2^' 

(3.22) ;=0 l=k+l 

1 _ log 2 

for some c, c independent of k (recall (|3.17|l l. Since the right-hand side in H3.22|l does 

not depend on k, we conclude that (|3.22|l holds for h G Ufc°=2['^fe"^ '^k'-il ~ '^i"'']- 
This completes the proof for = 1. 

In the case > 2 for a given (3 G (0, 1), let be the constructed above set 
and consider the direct product fi''. Let mes^ denote the Lebesgue measure in R'^. 
Recall that {fl'^)h = i^'^ — h, h G M.'^. Since for an arbitrary h = {hi, ■ ■ ■ , hd) we 
can write as a disjoint union = {0."^ n {il'^)h) U {Vl'^ \ {^'^)h) we have 

(3.23) mesd(r!'' \ [9.%) = mesd{fl'') - mes^f!'^ n {fl'')^). 

Noting that f^'^ n {n'^)h = (f^ n f^^ J x • • • x (f7 n flh^) and using IpO!^ for each of 
the factors we find 



mes^f^'^ n {n'')h) = Yl mesi(r! n n,,:^ 



(3.24) 

= n (mesi(r2) - mesi(Jl \ ^hj)^ 

By the construction of the set fi, there are ci, C2 > so that for \hj\ small, 

ci\hjf < mesi(f]\f],y) < C2\h^f 
which together with (|3.24|) implies 

d d 

(3.25) Yl (mesi(f7) - C2\hjf) < mesd(f7'* n {n'')h) < J| (mesi(r2) - ci\hjf). 

Recall that \h\ = ft. G M'* (we hope that denoting the Euclidean 

length of a vector and the absolute value of a number does not lead to confusion 
below). Note that 

d d 

(3.26) n (mesi(f^) - ci\h,f) = (mesi(f^))'' - ci ^ \h,f + E{hi, ■■■ ,hd) 

where the function E has the property that it is a finite sum of terms each of which 
contains a factor l^fei • • ■ |ftfcp |^ for some 1 < fci < • • • < /cp < d and p > 2. 
Note that each such factor is < 0[\h\^^) for small \h\. (Indeed, we can write 
\hi\l^\h2\'^ = . (|/ii|/|/i|)^ • i\h2\/\h\)'^ < \h\^'^ and for all other possible factors 
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again use the estimate \hj\ < \h\.) Substituting \A.'2(i^ in (|3.25|l and returning to 
(E2SIl, we find 

d d 

(3.27) ciJ2\hjf + 0{\h\^^)<mesa{n''\{n'')h)<C2Y,\h,f + 0{\h\^^) 

j=i j=i 

for small enough \h\ (we have used that mesd{^'^) — (mesi(ri))'*). We note finally 
that for any fixed /3 e (0, 1) there exist Ci, C2 > so that for all h 

d 

(3.28) Cil/il'^K^I/i.f <C2|/i|''. 

i=i 

(To see this, set h — h/\h\ and note that we just have to prove that for some 

Ci,C2 > 0, Ci < ^ ^2 if only J^Ui '^j'^ = ^- ^"^^^ 

the minimum and the maximum are attained at some points on the unit sphere.) 

Combining H3.27|l and (|3.28|) completes the proof of Lemma [2.41 
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